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Abstract 

We consider mirrors of the spherical shape, that can expand or contract. Due to the 
excitation of the vacuum around, some spherical waves radiated from vibrating mirrors are 
encountered. Using experience from well-known literature on studies of two-dimensional con- 
formal models, we adopt a similar framework to investigate such quantum phenomena in four 
dimensions. We calculate quantum averages of the energy-momentum tensor for s-wave ap- 
proximation. 
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1 Introduction 



According to Quantum Field Theory, a quantum vacuum state is the lowest eigenstate of the 
energy operator. In the language of the second quantization, it refers to the state with no 
particles. Sometimes, this simple basis state can become a quite complicated mathematical 
object if we try to describe it using observables related to quanta of elementary free fields. 
The fundamental observation is that all quantized fields have zero-point (vacuum) energies of 
size hu/2, where eigenfrequencies can be subject to some complicated equations defining the 
actual physical system. In other words, the physical vacuum can yield a complicated state 
involving virtual photons, electrons and other elementary quanta. Vacuum field fluctuations 
lead to significant physical effects, including charge and mass renormalization, Lamb shifts, 
Casimir effects or Unruh-Davies effects. The importance of vacuum fluctuations is recognized 
immediately whenever some non-trivial boundary conditions are assumed in field theoretical 
models. 

Among many interesting problems of this kind, some theoretical models are investigated where 
either a boundary of physical space or an interface between two different media changes its shape 
with time. A prototype example is a moving mirror (Unruh-Davies effect |ffl]). Induced by the 
motion of such objects, quanta of massless fields are created around by vacuum fluctuations 
and they can be radiated away. If a moving interface is neutral (as it usually occurs in physical 
applications), such radiation is weak. In practice, to get a significant radiation flux it is necessary 
to assume that interfaces or mirrors move with almost light-like velocities or undergo sudden 
or large accelerations. These theoretical ideas have been used to model the effects of strong 
gravitation on QFT or to derive Hawking radiation [Q, ||, to describe the squeezing effects in 
quantum optics 0], or to exploit the Schwinger suggestion how to explain the phenomenon of 
sonoluminescence Q. 

In this paper, we study the (3+l)-dimensional problem of a single spherically symmetric mir- 
ror. For the sake of simplicity, the case of massless scalar field is considered. The mirror acts 
as a perfectly reflecting infinite potential barrier, so this implies Dirichlet boundary conditions 
for quantum fields. Its radius depends on time, the mirror can expand or contract. We find the 
radiation flux. Our analysis is based on standard papers devoted to the two-dimensional moving 
mirrors || fj], |j, see also recent papers || [l(], 11, 12, 13, 14]. Maybe the most interesting 
physical situation involves a case when a mirror oscillates with the period related to the eigen- 
modes of the two-dimensional cavity. Under this resonance condition we get that the Casimir 
force is enhanced [15, [lfj. The main feature of two-dimensional situation is that the conformal 
symmetry allows one to transform the problem of time-dependent boundaries to the problem 
with static boundaries. It is not possible in the four-dimensional case. The quantum radiation 
of scalar particles by a moving plane mirror in the four dimensions has been considered in [17] 
(for non-relativistic motions). This paper addresses itself to the calculation of vacuum energy- 
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momentum tensor in the presence of a moving spherical mirror in four dimensions. However, 
we restrict ourselves to take account only spherical waves and spherically symmetric quantum 
excitations. 

2 Spherical mirror 

The model, we are dealing with in this paper, is the quantum theory of a massless real scalar field 
in (3+l)-dimensional Minkowski space-time. The field is almost free, no self-interactions are 
present except of interactions with a perfectly reflecting, spherical mirror. By a perfect mirror 
we mean a surface S on which the field is subjected to the Dirichlet boundary conditions, 

0| E = O. 

Center of the mirror rests in the chosen inertial reference frame while its radius R changes in 
time. 

As we want to compute the energy flux carried by the particles produced due to the motion 
of the mirror, it is crucial to have a well defined notions of the vacuum state and of the particle 
both in the past and in the future. It is possible if the mirror is initially and finally static, so we 
assume that its radius changes in time according to the formula 

{R for t < 
R(t) < t < h (1) 

Ri t>ti 

Due to spherical symmetry of the mirror, we are able to solve the d'Alambert equation, 

(-5, 2 + A)0(i,r) = O, (2) 
using the method of separation of the variables in the spherical coordinates. We get 

oo I 

cf>(t, r) = £ J2 C lm<t>l{t, r)Y lm (6, <p) + c.c. , (3) 

1=0 m=-l 

where Yi m (0,ip) are spherical harmonics, C[ m are arbitrary complex numbers and the functions 
4>i(t,r) are solutions of the equations 

d 2 1 d 2 1(1 + 

W2 --^r + A_Lij Mt ,r) = 0, Mt,W)) = 0, (4) 

with I being a positive integer. 

In the present paper we shall consider only the functions (j)(t, r) = 4>(t, r) which are inde- 
pendent of the angular variables. In other words, in the partial wave expansion (|3|) we confine 



2 



ourselves to the s-wave contribution with 1 = 0. This allows us to keep the technical details of 
the presented calculation at the relatively simple level. We plan to consider the more complicated 
case of the higher partial waves in the separate article. 

A general classical solution of the d'Alambert equation, assumed to be independent of the 
angles and to vanish on the surface of the mirror, can be written in the form 



<f>(t,r) 



du P^e-^fafar) + c.c, 



where 



Air I- > 



(5) 
(6) 



and (3{oj) is an arbitrary complex function. Ar and 5 R are defined by the formulae 

Ar = r - R , 5 R = 2 [R(t R ) - R ] . (7) 

R(t R ) is a position of the mirror in a retarded time t R , i.e. at the past moment when the wave 
incoming at the present time t and position r was reflected from the mirror, 



R(t R ) = t-t R . 



(8) 



Prom the eq. (0) it follows that t R and R(t R ) are functions of t — r. 

For t < the mirror is at rest and we can perform the standard quantization procedure by 
imposing on the field <j>i n (t,r) (now promoted to the level of an operator) and its conjugated 
momentum, 

Kin(t,r) = d t 4>in(t,r), 

the canonical, equal time commutation relations written in the spherical coordinates, 

i 



7r in (t,r),(f) in (t,r) 



;8(r 



0) 



Subscript "in" serves to recall that the relations above are valid only for t < 0. 

The field operator (pi n (t, r) can be expanded in the basis of classical solutions of the d'Alambert 
equation, 



4>in(t,r) 



1 



2ir Jo y/uj 



^(°)(r 



(10) 



where 



, fni . , sin (uAr) , , . 



for t < 0. 



The operators 6^,6^, satisfy the standard commutation relations (which follow from (0)), 

KM 



S(u — u'), 



(11) 
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The above relation allow us to interprete corresponding operators as an annihilation and creation 
operators for the field quanta. Consequently, we can define the state |0)j n which contains no 
particles through the equation 

b u \0) in = Vw. (12) 

From the moment t = 0, the mirror starts to move and finally, at t = ti, it reaches the shape 
of a sphere with the radius R = R\. In the period < t < t\ the evolution of the field operator 
differs from the free one in (^) due to destorsion of the modes caused by the interaction with 
the moving boundary, 

1 r°° clbJ r i 

cf>(t, r) = — -= b^Mt, r) + &t e **i/>Z(t, r) . (13) 



2ir Jn Ju 



While the system remains in the state \0)i n , it can no longer be viewed as containing no particles. 
This is due to the fact that for a generic motion of the mirror the functions ^(t, r) contain 
for t > an admixture of the negative frequency modes ~ e +lut , v > 0, and consequently 
cannot be interpreted as some annihilation operators for the "ouf vacuum. 

In order to calculate the amount of produced energy and momentum we shall compute the 
expectation value of the energy-momentum tensor T^ v . In the considered spherically symmetric 
case it has only two non-vanishing components: the energy density, 



2]\mj + {^ ■ (14) 



and the radial flux of the energy (or the momentum density) , 



_i(d<j>d^ d4>d<j>\ 

tr ~ 2 I &t 8r + drdtj • [ b) 

The corresponding continuity equation takes the form, 

d t T tt - -d r T tr = 0. (16) 
r 

Our main aim is to calculate quantum average values of the energy- momentum tensor. To 
establish a correct ultraviolet behaviour, we make use of point-splitting regularization method 

(6-0), 



(T t t) € = ^ in (0\d t <P(t,r)d t( j ) (t + ie,r)+d r <P(t,r)d r <P(t + i€,r)\0) m , 

(T tr ) e = ^ in (0\d t (l)(t,r)d r (l)(t + ie,r)+d r (l)(t,r)d t( l)(t + ie,r)\0) in . 

After straightforward calculations we obtain the following result for the energy density: 
ITV 1 ( 4 lr (XTOf v (l 2 \ A{l+TZ{u))[r-n{u)] H{u) \ 

\J-tt) = n — o I ~ o \U + ZK[U), UYq < 7. 77 „ , mo i > + 

x ' 32vr 2 r 2 I e 2 3 l w ' ss r\ e 2 + 4[r - TZ{u)} 2 1 + 2U(u) J 



(17) 
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(18) 



+ 4 Jin 



1 + 



A{r-lZ{u)f 



ln(l+27t(«)) 



and for the radial momentum density: 
1 



tr 



32vr 2 r 2 



^{u + 2n{u),u} s -^ 



1l(u) 



+ 



+ 0(e) 



K{u) 
1 + 2K(u) 



(19) 



where we have used an abbreviation lZ(u) = R(tfi(u)), u = t — r. The Schwartz derivative which 
appears in (^,19) is defined as 
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{u + 2K(u),u}, 



def 2 Tl (it) 



2ft(u) 



l + 27£(u) 2 \l + 2^(u) / 
As the formula (|i"8|) contains products of coinciding field operators taken in the same (for 



e — > 0) space-time point, it diverges in the limit of vanishing regulator. The appearance of 
this divergence reflects well-known ultraviolet problems in continuum quantum field theories. 
However, we should keep in mind here that the difference of energy densities between two 
physically realizable situations has a definite meaning and thus should be finite. 



Following this outline, we subtract from (18) the (regulated with the same prescription) 
vacuum energy density of the spherical wall with constant radius R\ , 



(Ttt, 



1 



Rl 32vr 2 r 2 I e 2 



! 1 . 
+ ^rln 



1 + 



4(r - -Ri 



2 4(r-Ri) 



r e 2 + 4(r - Ri) 2 



(20) 



and define 

{Tu) ren = 



lim((T tt ) 

e— »U \ 



(Ti 



ttIRi 



^(-\ {u + 2K{u)M S 



(21) 



2 J K(u)(r - Ri) + K{u) - i?i 
r 1 (r - R~ij(r - K{u)) 



+ 



+ 



In 



r - 7£(u) 



R 1 )(l + 2TZ(u)) 



1 + 272,(«) 

The expression (|l9|), which contains only products of distinct operators, is finite in the vanishing 
regulator limit, thus there is no need for any subtraction procedure, 



(Ttr)ren = ^MT tr Y ■ (22) 
e— >U 

The renormalized quantities (^j2^) satisfy the same continuity equation as their classical 
counterparts (14,15), 

dt(Ttt)ren d r (Tt r ) Ten = 0. (23) 

r 

It is worth to stress here that there is no contribution to the total energy from the subtracted 
expression (20). The quadratic term produces an infinite constant, while the other two terms in 
(pC|) return zero when integrated out over the whole space outside the sphere. 
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Using (21) we are finally at the position to calculate (for t > t\) the total energy. We obtain 
surprisingly compact formula, 

E = [ dV (T tt ) Tcn = — Heir ( ^ ) . (24) 

y t =const N ' 127T J R , \ 1 + 2K(U) ) V ' 



Let us note that respective boundary terms vanish here due to the assumption of static mirror 
asymptotical states and the symmetry of the considered physical system. We observe that the 
only term in Eq.([H|) which contributes to the total energy is just the one containing the Schwartz 
derivative. 

In the remainder of this paper, we consider an example being an explicitly solvable problem. 
It refers to the following time-dependence of the radius of the spherical mirror, 



R(t) 



R for t < 

t + R - a + aexp {-t/a) 0<t<i 

£ + JSo-&-*i + 6exp((ti-t)/6)) i<t<t x 

Ri t>t x 



(25) 



where 



U 



(Ri — R )t 



a(h - 1) 



(26) 



Here a, i, the initial and final radii Rq and R\ are the parameters defining the C 1 -class function 
R(t) and 

Vmax is a maximal velocity of the mirror 



1 — exp (—t/a). 



(27) 



The retarded time can be calculated by solving Eq.(H), 



tR(u) 



u + Rq for u < —Ro 
-a In [(a - u - Ro)/a] -Rq < u < av max - R 

h - bhx[(b + ti - u - R\)/b] av max - Rq < u < t\ - R\ 

u + R\ u> t\ — R\ 



The basic function TZ(u) is obtained immediately, 

Tl(u) = t R (u) - u, 
and allows to compute explicitly the total energy (p4|): 



E 



i?i — Rn 



247T a[R% - Rq + av max - i) 



l-v 2 



+ artanh v Tl 



(28) 



(29) 



(30) 
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3 Conclusions 



In this paper, the radiation emitted by the spherical mirror has been considered. We assumed 
the movement of the mirror lasts a finite period of time. This assumption helps to define and 
interprete asymptotic spaces of physical states uniquely. As the initial state it is considered the 
vacuum state (no particles present). The vacuum is perturbed by the moving mirror, and some 
flux of particles (radiation) is seen by an observer in the laboratory frame. The radiation is a 
pure quantum effect. We have calculated the vacuum expectation value of the energy-momentum 
tensor. To define finite results, we adopted the time-splitting regularization technique. To make 
calculations simpler, we restrict ourselves only to spherical waves and spherically symmetric 
quantum excitations. We obtained that the radiation flux depends only upon values of the 
mirror radius and its time derivatives evaluated along the intersection of the world history of 
the mirror with the observer's past light cone. 

4 Acknowledgements 

This work was supported in part by KBN grant No. 2 P03B 095 13. 

References 

[1] W.G. Unruh, Phys. Rev. D14 (1976), 870. 

[2] N.N.Birrel, P.C.Davies, Quantum Fields in Curved Spaces, Cambridge University Press 
(1982). 

[3] R.Brout, S.Massar, R.Parentani, PSpindel, Phys. Rep. 260 (1995), 329. 

[4] M.T. Jaekel, S. Reynaud, J. Phys. I France 2 (1992), 149. 

[5] J. Schwinger, Proc. Natl. Acad. Sci. USA 90 (1993), 958. 

[6] G.T.Moore, J. Math. Phys. 11 (1970), 2679. 

[7] S.A. Fulling, PCW. Davies, Proc. R. Soc. Lond. A. 348 (1976), 393. 

[8] P.C.W. Davies, S.A. Fulling, Proc. R. Soc. Lond. A. 356 (1977), 237. 

[9] V.P Frolov, E.M. Serebriany, J. Phys. A Math. Gen. 12 (1979), 2415. 

[10] M.Castagnino, R.Ferraro, Ann. of Phys. 154 (1984), 1; 161 (1985), 1. 

[11] W.R.Walker, Phys. Rev. D31 (1985), 767. 



7 



[12] R.D.Carlitz, R.S.Willey, Phys. Rev. D36 (1987), 2327; 2336. 

[13] G. Barton, C. Eberlein, Ann. of Phys. 227 (1993), 222. 

[14] G.Barton, Ann. of Phys. 245 (1996), 361. 

[15] C.K.Law, Phys. Rev. Lett. 73 (1994), 1931. 

[16] V.V.Dodonov, A.B.Klimov, D.E.Nikonov, J. Math. Phys. 34 (1993), 2742. 

[17] L.H.Ford, A.Vilenkin, Phys. Rev. D25 (1982), 2569. 



8 



